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Abstract 

Recently performed investigations of meson spectra allow us to determine the resonance struc- 
ture for the waves JJ^^ = 00++, 10++, 02++, 12++ = IJ'" = §0+ in the mass region up to 1900 
MeV, thus establishing the meson multiplets l^Pogg and 2^Poqq. Experimental data demonstrate 
that there are five scalar/isoscalar states in this mass region. Four of them are qq states, that is, 
members of the l^Poqq and 2^Poqq nonets, while the fifth state is an extra one not accomodated 
by qq systematics; it has the properties of the lightest scalar glueball. Analysis of the 00++-wave 
performed within the framework of the dispersion relation technique allows us to reconstruct the 
mixing of a pure gluonium with neighbouring scalar qq states belonging to I'^Poqq and 2^Poqq 
nonets: three scalar mesons share the gluonium state between each other - those are two compar- 
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atively narrow resonances /o(1300) and /o(1500) and a broad resonance /o(1530l25o)- The broad 
state is a descendant of the gluonium, keeping about 40-50% of its component. 

1 Introduction: retrospective view and the current state of 
the problem 

A great variety of the currently observed mesons and baryons represent systems built of quarks. 
These are baryons, which are three-quark systems (qqq), and mesons, which are quark-antiquark 
bound states (qq). More than 20 years ago the problem arose whether additional hadrons exist 
which are built out of another fundamental QCD particle, the gluon. An intensive search for the 
glueball - a particle consisting of giuons - was carried out throughout these decades. 

First evaluation of the glueball masses for different J^^ was done in the bag model Q . According 
to it, the lightest glueballs are scalars and tensors, 0++ and 2++; then follow pseudoscalar and 
pseudotensor glueballs, 0^+ and 2 

Recently considerable progress has been achieved in Lattice QCD calculations. The UKQCD 
collaboration obtained the following mass values for the lightest gluodynamical glueballs (i.e. the 
glueball without quark degrees of freedom taken into account): 

mG(0++) = 1549±53MeV, mG(2++) 2310±110MeV, mG(0"+) = 2332±264MeV. (1.1) 

Systematic errors are not included into the values given in Eq. (1.1); they are of the order of 100 
MeV. 

The IBM group obtained a slightly different value for the mass of the lightest scalar glueball [Q: 

"ig(0++) = 1740±71MeV, (1.2) 

The result of Rcf. g is as follows: 

tog(0++) = 1630±60±80MeV, mG (2++) = 2400 ± 10 ± 120 MeV. (1.3) 
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However, in Lattice calculations cited above the quark degrees of freedom have not been taken into 
account, since the existing computing facilities did not allow it. Quark degrees of freedom may shift 
the position of the glueball mass noticeably. The dispersion relation analysis of meson spectra jH , 
based on a restoration of the propagator matrix for scalar/isoscalar resonances, shows that the mixing 
with gg-states results in a mass shift of the order of 100-300 MeV. It should be stressed that, according 
to the 1/N expansion rules |^ [N — Nf — N^, where Nf and Nc stand for the light flavour and colour 
numbers), the mixing of a glueball with qq states is not suppressed. 

Experimental searches for the glueballs were particularly intensive during the last decade. There 
exist reactions where one could expect an enhanced production of glueballs. Central hadron production 
at the high-energy hadron-hadron collisions provides us with an example of such a reaction, because 
the particles in the central region are produced in the transition pomerons hadrons. The Pomeron 
is a gluon-rich system, so one could expect to find among secondary hadrons dominant production 
of glueballs, while the production of gg-states is expected to be small. However, the data on central 
hadron production, with statistics sufficient to perform a reliable partial wave analysis, are only now 
appearing. More accessible for experimental study appeared to be another reaction, which is governed 
by the transition gluons — > hadrons;this is, radiative decays J/i/' — > 7+hadrons . In these decays, 
hadrons are formed by gluons created in cc annihilation; therefore one may expect dominant production 
of glueballs in this reaction. Experimental study of hadron spectra in radiative J/^p decays has been 
carried out during 2 decades, and is still going on. Experimental information accumulated at the 
beginning of the 90's seemed to be rather discouraging, for in radiative J/ip decays qq states have 
been strongly produced: meson production branching ratios presented in the PDG compilation |^ 
show a number of resonances produced with similar probabilities, such as rj, rj' , /2(1270), /2(1525), 
etc., which certainly are gg-dominant systems. Such a situation presents a dilemma: 

(1) The glueball does not exist; it is "an unfulfilled promise of QCD" |lo| . 

(2) Glueball states are mixed strongly with the qq mesons, so in experiments one observes just 
these mixed states. 

Analysis of the 00"'"'" wave 0| definitely supports the second scenario. 

Experimental data on the transition form factors 77*((5^) — > Tr^,ri,T]' |pT| provide the following 
restrictions for probabilities to find the glueball components in rj and rj' mesons: Wrj < 8%, W^' < 20% 
This means that in the qq mesons observed in radiative J/ip decay one could find a glueball 
component at the level of 5-10%. Hence, the admixture of the qq component in the glueball should 
be considerably more, since the glueball can mix with several qq mesons. This qualitative estimation 
agrees with that obtained in the framework of the 1/iV expansion: according to this, the glueball 
component in each qq meson is of the order of 1/Nc, while the qq component in the glueball is of the 
order of Nf /Nc jl^. Of course, it should be stressed that some specific cases may differ from this 
general evaluation, because the mixing depends strongly on the relative spacing of mixed levels. 

If the scenario (2) with strong mixing of the glueball and qq states is realized in Nature, the 
search for the glueball is laborious and difficult work involving the identification of mesons and their 
systematics. Naive expectations, such as a study of the gluon-rich reactions with the purpose of seeing 
direct glueball production or, conversely, a study of hadron production in 77 — > gg-induced reactions 
with the hope of excluding gluon-rich hadrons, cannot be expected to succeed. 

The main channel of radiative J/?/' decays, as deduced from experimental data, is the production 
of broad hadron clusters. The production of these clusters may be viewed as a direct signal of strong 
mixing between the glueball and qq mesons. What happens is that, through mixing, one resonance 
accumulates the widths of other resonances. This effect has been first observed in where the 
low-energy part of the spectrum of the 00"'"+ wave has been analysed; this effect has been investigated 
in detail [t). When the two resonances mix completely with each other, one of them gets almost the 
whole width Fi + while the width of the other one tends to zero. In the case of an "ideal" mixing 
of three resonances, the width of one of them accumulates the widths of the other two, Fi + r2 + Fa, 
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and the widths of the others tend to zero. In reaUty, when the scalar gluebaU mixes with neighbouring 
states, there occurs a quahtatively similar effect; that is, a glueball situated among the scalar qq states 
mixes with them and accumulates a considerable part of their widths. From this point of view, the 
appearance of a broad resonance which is the glueball descendant is an inevitable consequence of the 
mixing. The broad resonance must be a neighbour of comparatively narrow resonances, which are 
the descendant of pure qq states; the broad resonance contains a considerable glueball admixture. 
The analysis of the 00++ wave in the mass range 1200-1800 MeV, based on the dispersion relation 
representation, reconstructs just this picture of the lightest scalar glueball mixing with qq members of 
multiplets I'^Pq and 2''Po- One may predict that such a scenario of mixing is common for all low-lying 
glueballs. 

Thus, the strong mixing of qq states with a gluonium does not allow easy identification of the 
glueball. In this case the only reasonable strategy is to study the systematics of all resonances in 
terms of qq multiplets. The extra states which do not fit into qq systematics should be regarded as 
candidates for the glueballs or other exotic mesons. This investigation program has been declared in 
[ p^ , and at the same time the first steps have been made in carrying it out: in Ref. iQ the /C-matrix 
analysis has been performed for the low-energy part of the wave 

IJPC = 00++. 

Detailed analysis of meson states in the region 1000-2000 MeV was possible due to the huge sample 
of experimental data collected in the latest decade by Crystal Barrel and GAMS Collaborations. The 
Crystal Barrel Collaboration has high-precision data on the production of three neutral mesons in 
the reaction pp annihilation at rest, 

pp {at rest) 7r°7r°7r° , 7r°7r°77 , 7r°?7?7 , (1.4) 

with the event numbers 720000 for (7r°7r°7r°), 280000 for {-K^TT^r]) and 185000 for (vr^ryry). The data 
on the reaction pp {at rest) —> Tr^Tr'^Tr", with somewhat lower statistics were published in 1991 [ p^ . 
However, the first fits of the spectra did not provide a correct identification of scalar resonances, for 
certain special features of the three-particle decay were not taken into account. A critical analysis 
of the situation has been made in 0, where it was shown that a resonance near 1500 MeV, 
which had earlier been identified as a tensor one, ^X2(1520), actually is a scalar resonance. Re- 
analysis of the reactions (1.4) within the T-matrix formalism performed together with the Crystal 
Barrel Collaboration has fixed the existence of new scalar resonances: /o(1500) ||l^ and ao(1450) 
[ p^ . In addition, in |l8| a noticeable production of the resonance /o(1360), with a half-width 

equal to 130 MeV, weas identified, although at that time it was not quite clear whether this was a 
newly observed resonance or a fragment of the broad resonance e(1300), which was widely discussed 
during the latest decades. Later on, after having performed the /^-matrix analysis for larger samples 
of data, it became clear that in this mass region there are two resonances - a comparatively narrow 
one, /o(1360), and a rather broad one , /o(1530l25o)- 

In the first stage of the investigation, the fitting to data was done in the framework of the T- 
matrix technique. The reason was obvious: the T-matrix representation of the amplitude is simpler 
for fitting to data; the advantages of the ii'-matrix approach reveal themselves only when there exists 
information on all possible channels of the reaction. In the mass region 1000-1500 MeV, there are the 
following channels in the 00++ wave: tttt, KK, rjrj and tttttttt, while in the region above 1500 MeV the 
channel rjrj' becomes important. It was obvious that the application of the sophisticated if-matrix 
technique for fitting to a limited number of channels (1.4) would lead to some ambiguities. 

The discovery of the resonance /o(1500) immediately gave rise to the hypotheses of its close relation 
to the lightest scalar glueball, and the possibility of such a relation was stressed in [|l^, In the 
following papers 1 13, 2^, |l[ |2^, several schemes have been suggested for the mixing of the lightest 
scalar glueball with the neighbouring gg-states. However, all these schemes did not take account of 
special features of the mixing which are due to the transition of a resonance into real mesons, although 
just these transitions, as was shown in a specified iC-matrix analysis, determine the structure of the 
00++-wave around 1500 MeV. 
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At the next stage of the 00++ wave analysis, the GAMS data on the spectra tt'^tt'^, 7777 and rjrj' 
have been included; these were obtained in the reactions |2^, |2^ : 

TT^p — > utt'^'k'^ , 717777 , 717777' , (1-5) 

together with the data of the CERN-Miinich collaboration : 

TT^P — > TITT+TT^ (1-6) 

and BNL |^ group: 

TTTT^Kk. (1.7) 

Simultaneous analysis of the whole data sample (1.4)-(1.7) was carried out in Refs. ||lj, ^ |3^, in 
the framework of the iiT-matrix technique; in this way the range of masses under investigation and 
the number of channels covered by the if -matrix fit of the 00+"'" amplitude gradually increased. 

The first investigation was done in the region of invariant meson energies y/s < 1100 MeV for 
the two channels only, namely, tttt and KK, In this analysis an observation, which became henceforth 
important, was made: the transitions which are responsible for the decay of meson states are also 
responsible for a strong mixing of these states. Moreover, the masses of mixed states differ essentially 
from the primary ones. These "primary mesons" were called in |l^ "bare mesons", in contrast 
to physical states, for which the cloud of real particles, tttt and KK, plays an important role in 
their formation. The masses of bare states are defined as the ii'-matrix poles. The above-mentioned 
accumulation of widths of the primary states by one of the resonances due to mixing, was also observed 
inRef. Q. 

As the next step, the if-matrix analysis was extended to 1550 MeV ||2^, with additional channels 
7777 and TTTTTTTr included. The channel Trmnr is rather important for the correct description of spectra 
from 1300 to 1600 MeV, since (T(7r7r TnnnT)/a{TnT tttt) is of the order of 0.5 at 1300 MeV and 



about 1.5 at 1500 MeV |31 . The use of channels tttt, KK and rjrj provides an opportunity to perform 
the qq classification of bare 00++ states, /o^'"''=, below 1600 MeV |||. The point is that qq-meson 
decays go to the new gg-pair via the production of intermediate gluons. According to the rules of the 
1/N expansion, the main contribution to the decay coupling constant comes from planar diagrams. 
When an isoscalar qq- meson disintegrates into two pseudoscalar mesons P1P2, namely, 

TTTT, KK, rjr] , rjri' , rj'jf , (1-8) 

the coupling constants can be determined, up to a common factor, by two factors. The first is the 
quark content of the qq-meson 

qq ^ nn cos (j) + ss sin(j) , (1-9) 

where nil — {uu + dd) / ^/2. The second is the parameter A, which characterises the relative probability 
to produce nonstrange and strange quarks by gluons in soft processes: 

uu: dd:ss ^1:1: \. (1.10) 

Experimental data provide the following values for this parameter: A ~ 0.5 js^ in central hadron 
production in hadron-hadron high-energy collisions, A = 0.8 ± 0.2 ||3^ for the decay of tensor mesons 
and A = 0.6 ± 0.1 ^ for the ratio of yields of 77 and 77' mesons in the decays J /ij) — > ^rj/^rj' . 

Coupling constants for the decay qq P1P2 into channels (1.8), which are defined by the leading 
planar diagrams in the 1/N expansion, may be presented as 

g(qq^FiP2) =Cp,p,(0,A)/, (1.11) 

where Cp^p,^{(j), A) is a wholly calculable coefficient depending on the mixing angle (j) and parameter 
A; is a common factor describing the unknown dynamics of the process. Therefore, experimental 
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investigation of resonance decays into channels (1.8) allows us to reconstruct the quark content of the 
state (i.e. its mixing angle (f)), thus making it possible to establish the meson systematics. 

However, on the basis of the decay constant analysis, it is impossible to determine unambigously 
whether we deal with a gg-meson or with the glueball. The reason is that the glueball decay is a 
two-stage process, with a subsequent production of two qq pairs. After the production of the first qq 
pair, in the intermediate state there exists a qq system with the following nn/ss content: 

nn cos 4>GluebaU + SS SlU (jjciueball , iB-'O- (j)GluebaU = \ ■ (1-12) 

For A — 0.45 — 0.80 the mixing angle is (j)Giuebaii =25° — 32°. At the second stage, the intermediate 
qq state (1.12) turns into the P1P2 mesons; this means that the relations between glueball coupling 
constants are the same as for the decay of the qq meson with = (jiciuebaii ■ 

Analysis of the tttt, KK and 7777 spectra performed in j29j proved that in the region below 1600 
MeV there are four scalar/isoscalar states, and only one of them is an ss-dominant state. Since each 
of the ^Poqq multiplets contains two J = states, which refer to two flavour combinations nn and ss, 
then, as a result of the analysis of the following dilemma becomes apparent: 

1. In the region 1000-1800 MeV there are three ^Poqq nonets: the basic one, l^Poqq, and two 
radial excitations, 2^Poqq and S^Poqq. In this case, there should exist two ss-dominant scalar mesons 
in the region 1600-1800 MeV. 

2. At 1600-1800 MeV, there is only one ss-dominant state. Then, one of the three mesons from 
the region 1200-1600 MeV is an extra one from the point of view of gg-systematics, and it should be 
considered as a candidate for an exotic meson: the ratios of couplings to the channels (1.8) found in 
[ p9| provide the basis to consider it as the lightest scalar glueball. 

Thus, after carrying out the analysis of Ref. ||2^, the immediate task was to extend the if-matrix 
analysis of the 00+"'" wave to the region 1600-1900 MeV. Such an extension suggested the inclusion of 
the 7777' channel into the fitting procedure; this has been done in [^0[ , where the /iT-matrix analysis has 
been performed in the mass region 500-1900 MeV, with the next five channels taken into consideration: 
TTTT, KK, rjr], tttttttt, rjrj' . It was shown that in the range 1600-1900 MeV there exists only one /o-meson 



with a dominant ss component, hence the analysis |30 confirmed the case 2. In this way, it was also 
shown that there are two variants for fixing the scalar glueball: 

Solution I. Two bare states, /i^""" (720 ± 100) and /o^^''<=(1260 ± 30), are members of the multiplet 
i^Poqq, and fl^''''%720) is the ss-rich state, with </)(720) = -69°±12°. The bare states /^^'■'=(1600±50) 
and /q '''''' (1810 ± 30) are members of the 2^Poqq nonet , and /o'^™(1600) is dominantly Tin-state, with 
</)(1600) = -6° ± 15°. The state /^'''■'=(1235 ± 50) is superfluous from the point of view of the gg 
classification; its coupling constants satisfy the ratios relevant to gluonium. Therefore, this state may 
be considered as a candidate for the lightest scalar glueball. 

Solution II. The basic scalar nonet is the same as in Solution I. The members of the next 
nonet, 2^Poqq, are as follows: /(^^■■<=(1235 ± 50) and /(^^"(1810 ± 30). Both these states contain a 
considerable admixture of the ss component: <^(1235) = 42° ± 10° and 0(1810) = -53° ± 10°. The 
state /g ^"(1560 ± 30) is an extra one from the point of view of gg systematics and it may be regarded 
as a good candidate for the lightest scalar glueball. 

The existence of two variants corresponds to the impossibility to answer unambiguously, on the 
basis of the information on the decay channels (1.8), if we deal with the glueball or gg-meson with the 
mixing angle (j> in the region 25° — 32°, as was stressed above. 

Both if-matrix solutions, I and II, lead to practically identical positions of the amplitude poles in 
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the complex mass plane. The amplitude has five poles: 

Resonance Pole position{in MeV) 



1015±15-i(43±8) 

1300±20-i(120±20) (1.13) 
1499±8-i(65±10) 
1750±30- i(125± 70) 
1530^250 -i(560± 140) . 

The broad resonance /o(1530l25o) ^^ot a new object in meson physics: this is that one which was 
called e(1300). A large width of /o(1530l25o) is due to the accumulation of widths of neighbouring 
resonances. 

Reliable and unambigous identification of the scalar glueball must be based upon the complete re- 
construction of the multiplets l^P^qq and 2^Paqq. Each of these nonets consists of two scalar/isoscalar 
states /o, one isovector/scalar state ao, and the scalar kaon Kq. As was stated above, it is reasonable 
to perform the nonet classification of highly-excited qq-states in terms of bare states, which do not 
contain clouds of real mesons. The analysis ||3^ fixed four /q"'"'^ mesons, which are necessary for 
the construction of two nonets; the two hghtest isotriplet resonances, ao(980) and ao(1450), are also 
known, see A full X-matrix analysis of the 10+"'" wave, |3^, provided the following resonance 
masses: 

ao(980) ^ (988 ± 6) - i(46 ± 10) MeV, ao(1450) ^ (1565 ± 30) - i(146 ± 20) MeV. (1.14) 

It should be pointed out that in the PDG compilation ^ the mass of the second resonance is too low 
by about 100 MeV. 

Corresponding bare states are as follows: 

a;j^'-<=(964± 16), a[;'^™(1670±70). (1.15) 

Identification of scalar resonances as members of the l^P^qq and ^^P^qq always raised problems. 
Namely, according to §,|3|l, the masses of the two lightest kaons are 1429±4±5 MeV and 1945±10±20 
MeV; these are noticeably higher than the average masses of other mesons which are candidates for 
the scalar-nonet members. This high position on the mass scale of the scalar kaon, ifo(1430), gave 
impetus to models where the basic I^Pq (j^-multiplet was fixed in the region 1350-1500 MeV, and the 
resonances /o(980) and ao(980) were considered as exotic states — hadron molecules [|8), multiquark 
bags HI], or minions 0, 

In Ref. [4^ the iC-matrix re-analysis of the S'-wave XTr-spectrum was carried out to determine 
the Kq^^'^. Another reason to re-analyse it was as follows. In Ref. the T^Tr-spectra have been 
investigated in two separate mass regions, 820-1580 MeV and 1780-2180 MeV, but the mass region 
1580-1780 MeV was not included into the analysis of the Ktt amplitude. Our experience in fitting the 
00"'"+-wave teaches us that separate consideration of different mass regions leads to the loss of 
certain information. In order to get a full picture, a simultaneous fit is needed; moreover, at 1580-1780 
MeV there is a rapid change of the amplitude. 

As follows from the ii'- matrix fit of the (IJ^ — ^0+) wave |Q, for a good description of the 
ii'Tr-spectrum in the region 800-2000 MeV at least two i^o-states are necessary. Correspondingly, the 
iO"'"-amplitude of this minimal solution has poles near the physical region on the 2nd sheet (under 
the i^TT-cut) and on the 3rd sheet (under Ktt- and i^T^'-cuts) , at the following complex masses: 

(1415±30) -i(165±25) MeV, (1820 ± 40) - i(125 ± 35) MeV. (1.16) 

The Kt]' threshold, being in the vicinity of the resonance (at 1458 MeV), strongly infiuences the 
iO"*" amplitude, so the lowest i^Q-state has a second pole which is located above the Kri'-cnt, at 



/o(980) 
/o(1300) 
/o(1500) 
/o(1750) 

/o(1530±«25o) 
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M = (1525 ± 125) - i(420 ± 80) MeV: the situation is analogous to that observed for the /o(980)- 
meson, which also has a two-pole structure of the amplitude due to the KK -Ihieshold. The K-q 
channel weakly influences the iO"*" Ktt amplitude: experimental data |]37| prove it as well as the quark 
combinatoric rules. 

The minimal solution contains two A'q^''° states: 

K'^^'^Umtfw) ' ATo "'=(18201^°) . (1.17) 

In this minimal solution, the lightest scalar bare kaon appears 200 MeV lower than the amplitude 
pole, and this latter circumstance makes it easier to build the basic scalar nonet, with masses in the 
range 900-1200 MeV. 

It is worth noting that the A'tt spectra also allow solutions with three poles, with much better 
still, for these solutions the lightest kaon state, Kq^'^°, does not leave the range 900-1200 MeV. In the 
tree-pole solution 

i^o^'^''"(1090±40), A:o'"'"(1375j:^25) ^ A'o^"'-^(1950t^°) , (1.18) 
and the i^Tr-amplitude has the poles: 

II sheet M = 998 ± 15 - i (80 ± 15) MeV 

II sheet M = 1426 ± 15 - i (182 ± 15) MeV 
III sheet M = 1468 ± 30 - i (309 ± 15) MeV 
III sheet M = 1815 ± 25 - i (130 ± 25) MeV. 

(1.19) 

The state A:^'"'^(1375j:if ), being near the Kij' threshold, results in doubling the amplitude poles 
around 1400 MeV. It should be underlined that masses of the lightest bare kaon states obtained by 
the two- and three-pole solutions coincide within the errors. 

The isT-matrix fit of the ^0+ wave makes it possible to complete, in terms of bare states, the 
construction of the two lowest scalar nonets. In line with the result for the 00++ wave |Q, where 
two solutions for an extra state (candidate for the glueball) were found, there are two variants for the 
nonet classification of scalar mesons. The basic l^Poqq— nonet is the same for both variants: 

a!j"''<=(960± 30) , /o^""(720± 100) , /o^"°(1260l Jo*') , K^^"" {1200^1%) . (1.20) 

It should be particularly stressed that the wave function /q'*''°(720) in the flavour space is close to 
the octet one; indeed, ^(720) = -69° ± 12°, while (/)octet = -54.7°. Correspondingly, /o'"^™(1260) 
is close to the flavour singlet. A similar situation is observed in the pseudoscalar sector where the 
flavour wave functions ij and rj' are close to the octet and singlet ones. It is even more analogous, if 
one takes into consideration that the mass diffrence of isoscalar states in these sectors coincide with 
each other, and the scalar masses are not much larger than corresponding masses of pseudoscalars, 
tTLscaiar — iTipseudoscaiar — (200±100) MeV. Such coincidcuccs clcarly point towards parity degeneration 
of the interaction forces in isoscalar channels. 

Thus, one may conclude: the basic nonet of scalar mesons is uniquely fixed by the _fC-matrix fit of 
meson sspectra. It is located rather low on the mass scale, in the range 750-1250 MeV. Here, at mass 
values below 1200 MeV, there is no room for exotic states. 

The 2^PQqq nonet contains the following states in Solution 1: 

a!j"'-^(1640±40) , /o^"°(1600±50) , f^""'' (ISlOtfoa) , 

K^''"'{1375tlf) or K^^'%1820tto) ■ (1-21) 
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An extra state with respect to the nonet classification is /q '^''''(1235l3Q°). 
In Solution II the 2'^Poqq nonet looks like: 

a^"'=(1640±40) , /o''"°(1235ii^") , /o''"<=(1810+?°o) , 

K}l'"'''{1375tlf) or K^""'" {1820tfo) ■ (1-22) 
In this solution an extra state is /g^'^°(1600 ± 50); once again it should be stressed that the mass of 
this state appears just in the mass region where Lattice calculations for the mass of the lightest scalar 
gluonium point; also the coupling constants with meson channels agree with the quark combinatorics 
ratios for the gluonium decay. 

Immediately after performing the X-matrix analysis in the range up to 1900 MeV, the problem 
of presentation of the 00++ amplitude as a dispersion integral has been raised. The dispersion N/D 
representation correctly restores analytic properties of partial amplitude in the whole complex s- 
plane. In addition, and this the principal point, within the dispersion representation it is possible to 
reconstruct the propagator matrix, thus evaluating the mixing of the qq-states and the glueball, and 
then to restore correctly the gluonium mass. The dispersive iV/I?-description of the wave 00++ has 
been performed in Refs. 0: in Rcf. the region 1200-1700 MeV where three scalar/isoscalar 
states are located, has been studied; then, in Rcf. Q the region under investigation has been extended 
to 1900 MeV, with the fourth state, /o(1780), included into consideration. 

The results of the A^/_D-representation of the 00++ wave allowed us to draw the picture of mixing 
for the lowest scalar gluonium: it mixes with the two neighbouring qq states — members of the 
multiplets I^Pq and 2'^Po: and the resonance- descendant of a pure glueball accumulates large parts 
of the widths of the neighbouring resonances, being transformed into the broad state /o (^^'^'^^^250) ■ 

It should be emphasized that the state fQ^'^°, which was found in the if-matrix fit, does not 
explicitly describe the gluodynamic glueball, for the state /g '^"^ contains nongluonic degrees of freedom 
related to real parts of the loop diagrams (imaginary parts are responsible for the decay process). The 
dispersion relation N/D method is able to restore the real and imaginary parts of the loop diagrams, 
thus providing a complete picture of mixing, so it also restores the mass of gluonium. In Solution I it 
is equal to 

rrigiuonium = 1225 MeV , (1-23) 

and in Solution II: 

mgiuonium = 1633 MeV . (1-24) 

The mass mgiuonium = 1633 MeV (Solution II) agrees well with the mass values of the lightest scalar 
glueball obtained in Lattice calculations. 

It is rather striking that both solutions obtained in the dispersive technique provide practically 
the same structure for the 00++ wave and the quark-gluon content of physical resonances: in both 

solutions the broad resonance, /o ^1530 +2®g*Jj^ , is a descendant of the gluonium, keeping about 40-50% 

of its component, while the remnant part of the gluonium is shared between /o(1300) and /o(1500). 
From this point of view, the structure of resonances in the range 1300-1600 MeV is uniquely solved. 

The formation of the broad state which is seen in the 00++ wave raises a question about the presence 
of such effects in other waves as well, for it is reasonable to believe that exotic mesons (glueballs and 
hybrids) with other quantum numbers can also afford a width accumulation. Because of that, the 
search for other exotic mesons must be inseparably linked with the study of broad resonances. 

There appears another problem which is not less intriguing: the broad resonance, after having ab- 
sorbed the widths of its neighbouring resonances, plays the role of a locking state, since it prevents the 
decay of neighbouring states with the same quantum numbers. This means that the broad resonance 
actually plays the role of a dynamic barrier for the nearby states. How does this dynamic barrier 
relate to the confinement barrier? — Only detailed investigations of broad resonances in other waves 
could answer this question. 
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2 iC- matrix and the dispersive A/'/D-representation of the scat- 
tering ampUtude 

In this Section a brief review is done for the technique used in the analysis of meson spectra. Namely, 
analytic properties of the amplitude are discussed together with the connection of the dispersion N/D 
representation to the K-matrix approach. The role of short and large distances in the formation of 
meson spectra under investigation is also discussed, and in this connection a notion "bare state" is 
introduced. The quark combinatoric relations between the decay couplings of the glueball to meson 
channels are considered in comparison with similar relations for gg-states: these relations provide the 
basis for the nonet systematics of mesons. 



2.1 Scattering amplitude, T- matrix and matrix 

Basing on a simple example, let us get over the terminology and notations used for the analysis of 

meson spectra. 

In terms of the wave function, which expresses relative movement of two spinless particles, the 
scattering at large distances is described by the incoming plane and outgoing spherical waves, with 
the coefficient f{6) related to the partial amplitudes as follows: 



oo 

f{e) = — ^(2^ + l)P,(cos e) [e^^^'^'^) - l] . (2.1) 



2ik 

This formula is written for the one-channel scattering in absence of absorbtion (fc is the momentum of 
relative movement, £ is the angular momentum and 6 is the scattering angle). The T-matrix element 
is determined by the scattering phase shift Sf -. 

Te = ^^[e^'^'-l]=e'^'smSi. (2.2) 
For the investigation of analytic properties, it is suitable to use an amplitude with other normalization: 

where p{k) is the invariant two-particle phase space factor: 

/I -73 1. 
dHP-M,k.), mP-MM = 2(^(^(2.)^^^^^(^-^-^^)- (2.4) 

Invariant phase space factor is determined by three four-momenta: the total momentum of scattcrred 
particles, P, with = s, and the momenta of particles 1 and 2, k\ and k^, respectively. For equal 
masses of particles 1 and 2, we get: 



p(.) = ^, k = ^^. (2.5) 
The ii'-matrix representation of the amplitude A(, reads: 

- l-ip{k)Ke{k-y 

is real in the physical region; the imaginary part of the amplitude is explicitly written in Eq. 
(2.6). In addition, Kg as a function of the invariant energy squared s is analytic near the threshold 
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singularity, s = Am?: a singular term is singled out, it is explicitly given by the two-particle phase 
space factor p. 

In presence of absorbtion, the scattering with an absorbtion is described by the absorbtion coeffi- 
cient rji inserted to the partial wave expansion (2.1): 

[e2»«^ - 1] ^ [r^^e^* - l] . (2.7) 

Here < ?7£ < 1; the case rj^ = corresponds to full absorbtion. 

It is suitable to display the energy-dependent amplitude Ti on the Argand-diagram, which is an 
appropriate instrument for searching resonances. The T-matrix element at fixed k (or s) corresponds 
to the point on the plane {ReT£,ImTi). As a function of k, it draws a trajectory on the circle with 
the radius equal to 1/2 and the centre in the point 0, i/2. In the inelasticity case, the trajectory 
enters the internal part of the circle. 

The if-matrix representation of the amplitude with an absorbtion requires fixing inelastic channels. 
Consider the inelasticity occurring due to another two-particle channel; we denote these channels by 
indices 1 and 2. Then the elastic scattering amplitude 1 + 1' ^ 1 + 1' (denoted as An, index i is 
omitted) can be presented in the form of Eq. (2.6): 

_ K{k^) 

1-zpiW)- ^^-^^ 
However, the block K(h?) has an imaginary part above the threshold of the second channel: 

K{k)^Kii+i- : — — . (2.9) 

1 - ip2J<22 

Here p2 is phase space factor of the second channel 2 + 2', and matrix elements Kn, K12 — K21 
and K22 are real functions of fc^ in the physical region. Threshold singularities of the channels 1 
and 2, which are located at s = (toi -|- m^)^ (threshold of the channel 1) and at s = (7712 -I- 1122)^ 
(threshold of the channel 2), are explicitly written in (2.8) and (2.9) ~ they are present in the phase 
space factors pi and p2, respectively. The function K{k'^) is real below the threshold of the second 
channel, (mi -I- m[)'^ < s < (m2 -I- m.2)^, because in this region p2 = i\p2\- 

2.2 Dispersion relation N/ D-method and the fC-matrix representation 

The dispersion relation N/D method [43j correctly reproduces analytic properties of the amplitude 
on the whole s-plane. Following ^|, we present here the elements of this method which are used 
in the analysis of meson spectra. 

Partial amplitude A{s) (as before, the index i is omitted for brevity) is written in form of the ratio 

^W = fg. (2.10) 

N{s) is a function of the complex variable s. It has only left-hand side singularities of the amplitude, 
which are related to the interaction forces, i.e. to the diagrams with meson exchanges in the crossing 
channels (see Fig. 1). These singularities are located on the left from the threshold singularities, at 
s = [mi + m[)^ - ml^^^^^„g. 

The iD-function contains only right-hand side singularities, which result from the rescattering of 
particles in the s-channel. Figure 2.1 shows the corresponding rescattering processes. 

First, we consider the one-channel case, with equal particle masses mi = m'l, in this case the 
D-function assumes the following form: 

r . ds' N(s')p(s') 
Dis)^l-Bis), B{s)= 7^^^- (2-11) 
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Here index 1 is omitted: toi — > m, pi — > p. The form of equation (2.11) suggests that D(s) — > 1 at 
s — !■ oo (more generally, D(s) — > Const at s — > oo, since this case can be reduced to Eq. (2.11) by 
the re-definition of N(s)). Moreover, in Eq. (2.11) it is also suggested that the D-function does not 
contain the Castillejo-Dalitz-Dyson poles (a detailed description of the iV/I?-method may be found 
inRefs. ||, ||). 

Representation of the iV-function in the form of a sum of the separable vertex functions is 
likely to be a reasonable anzatz: this technique is successfully applied for the description of the NN 
scattering amplitude ]45| ; in addition, the technique is elaborated for the presentation of the t-channel 
exchange diagrams as a sum of separable vertex functions p6| . For the simplest case, which is discussed 
below, N{s) = g^{s). Then 

Ais) ^ ^ . Bis) . r !^im^ . (2.12) 
1 — B[S) JAm^ Ti" s — s — iO 

Expanding Eq. (2.12) in a series with respect to B{s), we represent the amplitude A(s) as a sum of 
diagrams shown in Figs. 2a, 2b, 2c, etc.: B{s) in Eq. (2.12) is a loop diagram. At s > 4m^, the loop 
diagram is a complex quantity: 

ImB{S) = g\s)p{s) , ReB{s) T ^(^') 9^^s')p{s') ^^.la) 

J4,„2 TT S' - S 

The amplitude (2.12) stands for the case when the partial wave amplitude does not contain input 
particles: the bound states, if any, are formed by particle interaction taken in the iV-function. The 
inclusion of input particles into the amplitude corresponds to the assumption on D{s) increasing at 
s — !■ cx). The linearly growing D{s) can be written in the form: 

' - s — S — l\J 



4m^ 



The amplitude 



M^) = ^2. (2-15) 



^2 - s - B{s) 

is an infinite set of diagrams shown in Figs. 2d, 2c, 2f, etc.; B{s) stands for the loop diagram and 
(mg — s)~^ is the propagator of the input particle. 

The i^T-matrix representation of the A{s) amplitude is related to the explicit separation of the 
imaginary part of the loop diagram: 

^ K{s) ^ 9\^) (2.16) 

niQ — s — ReB{s) — ip{s)g'^{s) 1 — ip{s)K{s) ' — s — ReB{s) 

In the two-particle loop diagram, the function ReB{s) is analytical at the point s ~ Arn^. This means 
that poles are the only singulariries of K{s) in the physical region. However, at the left half-plane s, 
K{s) has singulariries related to the t-channel exchanges. 

The poles of the amplitude A{s), which are determined by the condition 

ml-s-B{s)=0, (2.17) 

are related to particles with quantum numbers of the partial wave under consideration. If the pole is 
above the threshold, at Re s > Am? , we deal with the resonance; this very case is studied further. Let 
the equality (2.17) be satisfied in the point: 

s^M'^ = p'^-iTp. (2.18) 
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Expanding the real part of the denominator (2.15) in a series near s — jj? , one has: 

ml-s- ReB{s) ~ (1 + ReB'{i/))(i/ - s) - ig^{s)p{s) . (2.19) 

The standard Breit-Wigner approximation comes when ImB{s) is fixed in the point s = fi^. If the 
pole is located not far from the threshold singularity s = Aw?, it is necessary to keep the s-dependence 
in the phase space factor, and we use the modified Breit-Wigner formula: 

^(^) = ^2_,1 i^p^,^ ' 7 = 1 + ReB\ij?) ■ ^^•^'^^ 

Similar resonance approximation may be also carried out for the /C-matrix amplitude representation, 
that corresponds to the expansion for K{s) given by Eq. (2.16) near the point s = /x^: 

K{s) = ^^ + f. (2.21) 



Here 

1 + ReB'ilj?) ' •' 2(1 + ReB'in'^)) 1 + ReB'{ii^) 



„2.^^_ i- g'(^') 2ff(/Z^)g'(M^) 



2.3 Multichannel scattering 

The resonance amplitude (2.15) can be easily generalized for the n channel case. The corresponding 
transition amplitude 6 ^ a is equal to: 

^<-i') ' J{-1'%) ■ 5(»)=|:S«M, (2.23) 

where B^c is defined by the standard expression (see Eq. (2.14)), with the properly chosen phase 
space factor, vertex function and integration region: 

g^{s')p{s') 9l{s')p,{s') , 4m2 ^ Ami (2-24) 

The transition amplitudes Aab form a matrix A. The amplitude written in the ii'-matrix representation 

reads: 

A = k ^ , (2.25) 

I-ipK 

where K is the n x n matrix, Kab{s) = Kba{s), I is the unit nx n matrix I = diag{l, 1, . . . , 1), and p 
is the diagonal matrix of phase space factors: 

p = diag{pi{s), p2{s), p„(s)) . (2.26) 

The ii'-matrix elements are equal to: 

.9a(-s).9b(g) „„^ 

^''"(^^ = ml-s- ReB{s) ' ^'-'^^ 

In the vicinity of the resonance, the JC-matrix elements may be expanded in a series: in this case we 
have a representation of the iC-matrix elements similar to that of Eq. (2.21). 
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2.4 gg-mesons: the problem of small and large distances 

The qq classification of meson states in ttie vicinity of 1000-2000 MeV faces the problem of the 
quark-hadron duality as well as a tightly related problem of separating large- and small-r interactions 
contributing into formation of meson spectra. 

Let us discuss these problems using the language of the standard quark model. In this model 
the qq levels are determined by the potential which increases infinitely with r: V{r) ^ ar (see Fig. 
3a). Infinitely rising potential creates an infinite set of the qq levels. However, it is obvious that the 
standard quark model is a simplified picture, since only the lowest qq-\eve\s are stable with respect to 
hadron decays. The heavier states are decaying by hadron channels: an excited {qq)a-state produces 
new qq pair, then the {qq)a + (qq) quarks recombine into mesons, which leave the confinement trap 
with a formation of the continuous meson spectrum. This structure is conventionally shown in Fig. 
3b, where the interaction related to the confinement is shown as a certain potential barrier, namely, 
the interaction at r < Rconfinement creates the discrete levels of qq spectra, while the transitions 
into the region r > Rconfinement givc rise to the continuous meson spectrum. This very meson 
spectrum is observed in the experiment. This problem of reconstruction of the gg-levels created at 
r < Rconfinement IS dircctly related to the determination of the effect of the meson decay spectra on 
the level shift: classification of # levels requires elimination of the components of real mesons that 
are the decay products of real meson. 

The iiT-matrix representation of the amplitude resolves the problem of excluding the components 
of real mesons; formally, it corresponds to the limit — > in Eq. (2.25). If only leading pole 
singularities are taken into account, the transition amplitude 6 — s- a assumes the form: 



Thus, the if-matrix pole corresponds to the state with the removed cloud of real mesons. For this 
reason, we call corresponding states "bare mesons" ||2^, However, one should distinguish between 
this notation and that of " bare particles" of field theory, where a cloud includes virtual mass-off-shell 
particles as well. 

In the case when the gg-spectrum contains several states with the same quantum numbers, the 
amplitude ^l^°[s) is determined by the sum of the corresponging poles: 



Representation of the amplitude responsible for the interaction at r < Rconfinement in the form 
of a series of poles is not new: it was widely used in dual models for contributions which are leading 
in the 1/iVc-expansion. From the point of view of these models, the s-independent term fab is the 
thorough contribution of poles which are distant from the region under consideration. 

Coupling constants of the bare state, gi'^^K), are a source of information on the quark-gluon 
content of this state. 

2.5 Coupling constants of the gluonium and gg-states to meson channels: 
the rules of the 1/A^-expansion and quark combinatoric relations 

Quark-gluon content of states related to the if-matrix poles (bare states) is revealed in relations 
between couplings of these poles to meson channels, gi°'\ 

First, let us evaluate these couplings using the rules of 1/iV-expansion; this evaluation will be done 
both for the transitions glueball — > two mesons and for the transitions qg— state two mesons. 



ga{K)gt{K) 
fi^ — s 



+ fab ■ 



(2.28) 




(2.29) 



13 



For this purpose, we consider the ghion loop diagram which corresponds to the two-gluon self-energy 
part: glueball two gluons glueball (see Fig. 2.4a). This loop diagram is of the order of unity, 
provided the glueball is a two-gluon composite particle: B{G gg ^ G) ^ 9G-*gg-^c ^ where 
gc^gg is a coupling constant of a glueball to two gluons. Therefore, 

ga^gg ~ l/N, . (2.30) 

Coupling constant for the transition gc^qq is determined by the diagrams of Fig. 4b. Similar evalu- 
ation gives: 

gG~.qq ^ gC^gggQCD^c ~ 1/iVc . (2.31) 

Here ggcD is the quark-gluon coupling constant, which is of the order oil/^/Nc Coupling constant 
for the transition glueball two mesons in the leading 1 /Nc terms is governed by the diagrams of 
Fig. 4c: 

gh^mm - 9G^gggrn^,qNc ^ 1/iVc • (2.32) 

Here the approximate equality gm—>qq ~ 1/ \/N~c is used which follows from the fact that the loop 
diagram of the meson propagator (see Fig. 4d) is of the order of the unity: B{m qq ~^ m) 
9m->qq^c ~ 1 • The diagram of Fig. 2.4e-type governs the couplings for the transition glueball 
two mesons in the next-to- leading terms of the 1/A^c-expansion: 

gatmm ^ 9G^gggQCD9m^qqNc ^ 1/^c • (2-33) 

As was mentioned above, the glueball can decay into channels (1.8). In this way, the production 
of light quarks by gluons is going on with a violation of the flavour blindness, uu : dd : ss — 1 : 1 : X. 
Within such an assumption, the coupling constants for the transition glueball two pseudoscalar 
mesons can be calculated using quark combinatorics. The rules of quark combinatorics were suc- 
cessfully applied to the calculation of yields of secondary particles in the hadron-hadron high energy 
collisions |Q and to the decay J/ip ^hadrons [Esl. The calculation of couplings in the decays 
glueball mesons has been carried out in Refs. |20| , |3^ , 

The glueball decay couplings to channels (1.8) are given in Table I for leading, g^^^^, and next- 
to-leading, gQ^jnrn, terms of the 1/iV-expansion. Unknown dynamics of the decay is hidden in the 
parameters Gl and Gnl- The decay constant to the channel n is a sum of both contributions: 

9h^mM + gatrnM . (2.34) 

The second term is suppressed, as compared to the first one, by the factor N^, the experience in 
calculation of quark diagrams teaches us that this suppression is of the order of 1/10. 
The sum of couplings squared satisfies the sum rules: 

E„ {9h^n.Mfl{n) = \Gl{2 + Xf , E„ {9^trnJn)f lin) = iG^ J2 + A)^ , (2.35) 

where I{n) is an identity factor for the particles produced, see Table I. These sum rules follow from 
the quark-hadron duality: the sum of squared coupling constants over the whole set of flavour states 
is equivalent to the sum of cut diagrams with the quark loops (diagrams of the type of Fig. 4f for 
leading and of the type of Fig. 4g for next-to- leading terms). Each quark loop contains the factor 
(2 -f A) related to the summation over light flavours, see Eq. (1.10). 

Quark combinatoric rules may be applied to the calculation of couplings of (<7g)a-mesons to pseu- 
doscalar channels (1.8). There exist two types of transitions (55)0— state^two mesons: they are 
shown in Fig. 5a and 2.5b. The type of process represented by the diagram of Fig. 5a is leading 
according to the rules of 1/iV expansion; its coupling constant is of the order of 

9m{a)->mm 9m^qq^c ^ ^= ■ (2.36) 
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The decay constant for the process of Fig. 5b is of the order of 



„NL ^ 3 A^2 4 ^ ^ (r, 

The couphng constants for the decays {qq)a tttt, KK, rj-q, rji]' , rj'rj' are shown in Table 2 both for 
the leading and next-to-leading orders; and g^^ are parameters, in which the unknown dynamics 
of soft decay is hidden. Concerning the glueball decay, the coupling for the {qq)a-meson decay to the 
channel n is the sum of both terms: 

9m(a)^mmi^) + 9m{a)^mmi^) ■ (2.38) 

These two terms in Eq. (2.38) define the decay constant of the (gg)a-meson in general case: different 
variants of fixing the ratios of coupling constants correspond to different choices of the ratio Ql/ 9nl- 
The examples of different fixations of gh/gNL may be found in p9[ . 

Let us stress once again that the coupling constant ratios for the {qq)a-st&tcs (Table II) become 

A. 



identical to those of the glueball, when (j) = arctany-^: it is valid for leading and next-to- leading 
contributions. Therefore, that, on the basis of a study of couplings to the hadron decay channels only, 
it is impossible to distinguish between the glueball and the I — Q (gg)a-meson with the mixing angle 
6 close to 30°. 
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Fig. 1. Complex s-plane and positions of singularities of the partial amplitude: right-hand singu- 
larities at s > {mi + m'lY' are due to the elastie and inelastic rcscattcrings, left-hand ones are due to 
the interaction forces, that is, particle exchanges in the crossing channels. 
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Fig. 2. Diagrams standing for the s-channel scattering. 
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Fig. 3. a) Standard quark model potential with stable qq-\cvc\s, b) Potential with the unstable 
highly excited levels which corresponds to the realistic situation for qq states. 
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a) b) 
Fig. 5. Diagrams for the decay of the {qq)a-sta,te into two mesons. 
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3 A'-matrix analysis of meson spectra and nonet classification 
of gg- states 

In this Section, the resuhs obtained in Ref. Q |36[ || for the waves 00++, 10++ 02++, 12++ and 
|0+ are presented. On the basis of this analysis, the nonet classification of qq-states is established. 

3.1 i^'-matrix fit of 00++-wave: the spectra tttt, KK, rjr] and r]r]' 

To describe the spectra, in Ref. |^ the standard if -matrix representation of the 00++ amplitude 
(2.25) was used, where Kab is the 5x5 matrix (a, 6 = 1,2,..., 5), with the following notations for 
channels: 1 = tttt, 2 = KK, 3 = 7777, 4 — rjrj' , 5— {71711171+ other multimeson states). 
Matrix elements of Kab are parametrised in the form similar to Eq. (2.28): 

V a ^o^-s s + so J 

with the restriction sq > 1 GeV^. 

The following formulae provide the description of the spectra tttt, rjrj and 7777', obtained by GAMS 
for the reactions with the i-channel pion exchange: 

A^N^Nb - gi^Nl5^N)F{t)D{t) Kia{t)il - ipK)ab\ (3-2) 

^f..^ + /,.(0i^^). (3.3, 

y - s s + so J 

Here D{t) is the pion propagator, F/v(i) is the nucleon form factor related to the vertex ttNN, and 
g['^\t) and fia{t) are form factors of the pion block. 

The piece of the amplitude for the reaction pp {at rest) — s- 7r°7r°7r°, 7T°riri, which describes the 
production of two mesons in the 00++ state, is written as follows: 

App — *mcsons = Ai{s23) + A2{si3) + Asisu) (3.4) 

The amplitude Ak{sij) corresponds to the process with the "last interaction" of particles ij, while the 
particle k remains a spectator. 

The amplitude ^1(323) for the spectra tttt and 7777 has the following form (6 — tttt, 7777): 

^1(S23) = Kpp^,a{s23) (1 ^ ipK)ab , Kpp^^a{sij) = ( X! A/fr^" + 'I'PP^-'^ 1 GcV + Sq \ ^ ^^^^^ 

In the reaction pp {at rest) — > 7r°7r°7r°, the amplitude is symmetric with respect to the permutation of 
pion indices: Ai{sij) = ^2(5^) — ^3(sy ). The 7r°7r° interaction in the reaction {at rest) 7T°7i°ri 
is determined as follows (below 6 = tttt): 

Al{S23) ^ Kpp,^.a{S23){l~ipK)^^ , Kpp^^a{Sij) = \ JTYZ + <Pppri,a J ■ 

\ a " " / 

(3.6) 

Parameters A^^^^, 0pp7r, ^pp-q and tppprj can be complex quantities with different phases, that follows 
from the three-particle interaction; more detailed discussion of the amplitude pp ~^ three mesons is 
given in Ref. ||50|. 



19 



3.2 Results of the iT- matrix fit for the 00 
MeV 



wave in the region below 1900 



'9k 



Simultaneous if-matrix fit of tlie 00++ spectra in the mass region 550-1900 McV performed in 
pointed out to the existence of five bare states, /o^"°. Only two of them, /|5'^''=(720) and /o''"'=(1810), 
contain a large ss-component. This means that only two '^Po'Z^-nonets can be built in the mass region 
below 1900 MeV. 

The following requirements provide the ground for the nonet classification of bare states: 

(1) Nonet partners are orthogonal in the flavour space, i.e. they must have the mixing angle differences 
(see Eq. (1.9)) equal to 90°: - 02 = 90° (for this value the corridor 90° ± 5° is allowed). 

(2) Coupling constants and g^^ (see Table 2) are approximately equal for nonet partners: 
and g^^ g^^. 

The standard quark model requires the equality of coupling constants. However, the s-dependent 
vertex functions and loop diagrams violate this equality because of the presence of mass differences 
for nonet partners. Moreover, the ii'-matrix coupling constants have an additional s-dependent factor 
(1 + B'{s))~^ (see Eq. (2.22)): this factor strongly affects the region of small masses (the region of 
the basic nonet l^Po)j where the thresholds and left-hand singularities of partial amplitude play more 
significant role. 

The fitting to experimental data (1.4)-(1.7) resulted in two solutions, I and II. First, let us sum 
up the results for Solution I: 



Type of the state 



jbare(720) 
ybare(j260) 

/□^'''■^(leoo) 

/iJ'*'-'=(1810) 
/o^"''^(1235) 



2^Poqq 
Glueball 



Flavour wave function 
0.40nn- 0.92ss 
0.92nn -h 0.40SS 
0.995nfi- O.lOss 
0.10nn-H 0.995ss 
0.91nn -I- 0.42ss. 



(3.7) 



In Eq. (3.7) the "flavour wave function" is introduced for the glueball: it describes the flavour content 
of intermediate state for the glueball decay: see Fig. 4c. 
Now let us summarize the results for Solution II: 



Type of the state 


Flavour 


wave function 


y bare (720) - 


^ l^PoW 


0.40nn - 


- 0.92SS 




l^PoW 


0.92nn- 


1- 0.40SS 


/o^'^'-'=(1235) - 


23PoW 


0.74ml - 


1- 0.67SS 


/o^'^'-'=(1810) - 


23PoW 


0.67nn - 


- 0.74SS 


/o^'^'''=(1600) 


Glueball — s 


■ 0.91nn- 


1- 0.42SS. 



(3.8) 



The quality of data description by Solutions I and II can be seen in Figs, 
curves, respectively). 



6-3.4 (dashed and solid 



3.3 The resonances: are they bumps or dips in the spectra? 

During decades the search for meson resonances meant the search for bumps in the particle spectra. 
Only recently the understanding came that it is not always so, and the resonance /o(980) provided 
us with an example. In the peripheric vrTr-spectra it reveals itself as a dip, and a set of papers was 
devoted to this phenomenon, see |9[ 0. The study of 00++ wave proved that meson resonances in 
the region 1000-1600 MeV appear not only as bumps or dips, but as shoulders in the spectra as well. 
The main particular feature of the resonance is not a bump or a dip in the spectrum, but a circle on 
the Argand diagram. 
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Figures 10 and 11 demonstrate the Argand diagrams relevant to the fits of spectra under discussion. 
In Fig. 10 one can see the 00'*"''" amplitudes A-^tt^tttt, A^-k^kki A-mr^riri and ^tttt-^jjjj' as functions 
of the energy. Indeed, there exist rather legible circles obtained for the resonances /o(980), /o(1300), 
/o(1500) and /o(1780). The manifestation of resonances /o(980) and /o(1300) in the form of circles is 
rather clear for the amplitudes at large \t\, see Fig. 11. 

3.4 Resonance /o(980): is it the KK-molecule or the descendant of the 
Hghtest scalar gg-states? 

This is a principal problem for the qg-systematics that was firstly investigated in by using the K- 
matrix analysis for the low-energy part of the 00"'"^-wave. Following this paper, we provide arguments 
that /o(980) is the descendant of the lightest qq-states. 

Two poles correspond to the resonance /o(980), with the following complex masses: M ~ 1015— i46 
MeV on the second sheet (under the tttt cut) and M — 936 — i238 MeV on the third sheet (under the 
TTTT and KK cuts). The second pole, at 936 — i238 MeV, appears because of the well-known doubling 
of poles affected by the proximity of the ii'A'-threshold - for example, see Refs. [|l|, |2|. The first pole 
at 1015 — i46 MeV dominates in the TTTr-spectrum, providing a sharp dip in the tttt ^ tttt spectrum 
or a bump in the spectrum Tnr{t) — s- tttt at large \t\ p^ . Let us study the dynamics of this pole when 
the decays into the channels 1 = tttt and 2 = KK are subsequently switched on and off. To this end, 
let us make the following replacement in the i^-niatrix amplitude: 

5i(720)^affi(720), .92(720) -> 652(720), (3.9) 

with the parameters 6 varying in the interval < 6 ^ 1 ■ At 6 ~* and ^2^0 the decay channels 
for the lightest OO^"*" state are switched off, and we deal with the "bare" state, in terms of Refs. 
[|9[ At 6 = '^2 = 1 the real case is restored. At 6 — 6 — the mass of the bare state is in the 
vicinity of 720 MeV, while, with the increase of ^a, the lightest scalar state acquires the components 
of real mesons, tttt and KK, and, due to transitions to these states, it mixes with the other scalar 
ones. As a result, the mass of the lightest scalar state increases, approaching the region around 1000 
MeV. At — ^2 — ^ the pole of the amplitude is at 

M{real position) = 1015 - i46 MeV , (3.10) 

i.e. near the iiT ^-threshold. Therefore, the iiT .R'-component is of the quasi-molecular type: the relative 
momentum of iiT-masons is small, hence the mean distance between the if-mesons is comparatively 
large. However, one can see that the if.R'-component weakly affects the formation of the final state. 
Indeed, let us switch off the KK-sl&ie, i.e. let us put in the amplitude 6 = 1 and 6 = 0. Then, the 
pole appears in the point: 

Miwithout KK) = 979 - ibi MeV . (3.11) 

The mass shift 

M {real position) - M {without KK) = 36 + i7 MeV (3.12) 

is actually not large, thus proving the role of the ii'A'-component in the formation of the real /o(980) 
state be insignificant. 

Concluding, the 7^-matrix analysis of 00+"'"-wave restores the following picture of the formation of 
/o(980). Before mixing, there existed the lightest scalar/isoscalar gg-state /q°'''^(720± 100), with the 
flavour wave function close to the octet one: 

V'/ia™nr(720) = CosOg ^8 " SmO s ^1 , 0s = 14° ± 12° , (3.13) 

1 - 1 - [2 

ijji = —^{uu + dd + ss) ips = —i={'uu + dd) — \ -ss . (3.14) 
V3 v6 V 3 
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The mixing with other states, which is going through the transition /q'^'''^(720) — > tttt, leads to the 
formation of the resonance with the characteristics which are almost the same as those observed at 
the experiment, see Eq. (3.11). The onset of the _ft'A'-component /q'"''^(720) KK affects relatively 
small shift of mass and width, see Eq. (3.12). 

It should be noted that direct measurements also point out that /o(980) has considerable small-r 
components: the production of /o(980) is not suppressed in the reaction n^p /o(980) p at large 
momentum transfers ^ as well as in radiative J / ip decay |Q . 

One should pay attention to the fact that the lightest scalar/isoscalar state, that is, 77- meson, also 
has flavour wave function close to the octet one: 77 — CosOp i/jg — SinOp ipi with 9p = —16.7° ± 2.8° 

[1- 

3.5 The wave IJ^^ = 10++ 

Two isovector/scalar resonances are clearly seen in the reactions (1.4) [|8[ |l9|. The lightest of them 
is a well-known ao(980), while the next resonance is the recently discovered ao(1450): according to 
PDG compilation §, its mass is 1450 ± 40 MeV and the width is T = 270 ± 40 MeV. Note that the 
fitting to the latest high statistics data ^ |5^ provided us with a greater mass value — it is equal 
to 1520 ± 40 MeV. 

For the description of the scalar/isoscalar amplitude, in Rcf. the two-pole 4x4 matrix 

was used, with the channel notations 1 = ttt], 2 = KK, 3 = Tryy', 4 = multiparticle states. The 
couplings to the two-meson channels are defined by the quark combinatoric relations, see Table 3: we 
remind that the constants are the same for all multiplet members. At the first stage of the fit, the 
coupling constants of the lightest resonance oq have been varied in the interval limited by the constants 
^Ljj6are(720)] and 5^ [/o*°'''= (1260)] . In all the variants of the fit, the coupling g^[a^„'''''' {lightest state)] 
was found to be rather close to g^[fQ°'^'^{126Q)]: in the final variant of the fit these coupling constants 
are taken equal to each other. The couplings of the next isovector/scalar resonance to two mesons are 
also fixed: they are equal to each other for all the multiplet members 2^Po. 

The fit allowed us to find two solutions for the wave 10++, which practically coincide for the 
resonance sector and differ for background terms. Positions of the amplitude poles and relevant bare 
states are shown in Eqs. (1.14) and (1.15). 

In Ref. the hypothesis was discussed about resonances ao(980) and /o(980) belonging to a 

special class of states (minions), which are loosely bound to hadron channels: the small widths of 
ao(980) and /o(980) were considered as arguments in favour of that particular nature. The character- 
istics of ao(980) provide a good opportunity to check this hypothesis, because actually the components 
of real mesons yielded by the decays ao(980) nrj,KK do not influence this state, see Eqs. (1.14) 
and (1.15). However, the if-matrix data fittings proved that g'^[/o"'^(964)] is not small, being of 
the order of a standard hadronic value, for the small width of ao(980) is related not to the probability 
of the decay, as it could follow from the minion nature of ao(980), but to the threshold effect. It 
should be emphasized that this result is also seen in the T-matrix analysis of the data M, M . 



22 



Tabic 1 

Coupling constants of the glueball decaying to two pscudoscalar mesons, in the leading 
and next-to leading terms of 1 /N expansion. O is the mixing angle for Tj — rf mesons: 

rj = nn cos © — .s.s sin and r/' = nn sin Q + .s.s cos 0. 



Channel 


Glueball decay 
couplings in the 

leading term of 
expansion. 


Glueball decay 
couplings in the 

next-to-leading term 
of 1 /A^ (expansion. 


Iden- 
tity 

factor 


n n 







1/2 


7r~'~7r 




n 
u 


1 

± 


K+K- 







1 


K°K° 


VXG^ 





1 


rjrj 


G^ (cos^e + Asin^ e) 


2G^^(cose-yjsine)2 


1/2 


m' 


G^(l- A)sine cose 


2G^-^(cose- y^sine)x 


1 






(sine + y^cose) 




in' 


G^ (sin^ e + A cos2 e) 


2G^^ ^sine + y^cose^ 


1/2 
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Table 2 

Coupling constants of qq-meson decaying to two pseudoscalar mesons in the leading and 
next-to-leading terms of the expansion, (j) is the mixing angle for nh and ss states, see (1. 



Channel 


Decay couplings of 
qq mesons 
in leading term 
of 1/N expansion. 


Decay couplings of 
qq mesons 
in next-to-leading term 
ofl/N expansion. 




cos(j)/\/2 







TT+TT" 


cos0/a/2 







K+K- 


5f^(\/2 sin \/A cos ?i)/\/8 









5^(V2sin(?!) + yXcos (j))/V8 







m 


g^ (cos^ e cos(^/\/2+ 

\/A sin (j) sin^ 0^ 


V2g'^^ (cos e - sin 6) X 
(cos ^ COS 9 — sin sin 9) 


m' 


g^ sin 6 cos Q (cos $/a/2- 


V 2^ 


(cos 9 - -y^ sin9)x 






(cos sin 9 -|- sin (p cos 9) 

+ (sin9 + cos9)x 
(cos ^ sin 9 — sin (j) cos 9)] 


ri'ri' 


g^ (sin^ e cos (j)/V2+ 

\/A sin $ cos^ 9^ 


v^g-""'^ (sin 9 + cos 9) x 
(cos COS 9 -|- sin (/) sin 9) 
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Fig. 7. Evc;iit numbers versus invariant mass of the TTTT-system selected for various intervals of the 
momentum transfer squared t. Solid curve stands for Solution II, dashed one to Solution I. 
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Fig. 8. Description of the angle moments for the tttt distributions measm'ed in the reaction 
TT-p 7111-71+ [p7|; Solution II |a. 
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Fig. 9. The 7r°7r° spectra in the reactions pp Tr^Tr^Tr", pp — > ryTr^Tr" and rjr] spectrum in the 
reaction pp —> ir'^rjr] and rjir^ spectrum in the reaction pp — > tt^tt^?]. Curves correspond to Solution II 

m. 
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Fig. 11. Argand diagram for the isoscalar S-wave 7r7r(f) mr scattering amplitude |30 at different 
momentum transfers squared, t. 



3.6 K-matrix analysis for the Ktv S-wave 

Partial-wave analysis of the K^w^ system for the reaction K^p K^-rr^n at 11 GeV/c was carried 
out in Ref. where two alternative solutions (A and B), which differ only in the region above 1800 
MeV, were found for the 5- wave. In that paper the T-matrix fit for the Kt: S-wave was performed as 
well, though independently for the regions 850 — 1600 MeV and 1800— 2100 MeV. In the first region, 
the resonance -K'o(1430) was found: 

Af/j = 1429 ± 9 MeV, r = 287 ±31 MeV. (3.15) 

In the second mass region, Solutions A and B provided the following parameters for the description 
of the resonance ^^0(1950): 

Solution A Mr = 1934 ± 28 MeV, T = 174 ± 98 MeV, 

Solution B Mi^ = 1955 ± 18 MeV, T = 228 ± 56 MeV. (3.16) 

The necessity to improve this analysis is obvious. First, the mass region 1600 — 1800 MeV, where 
the amplitude varies quickly, must be included into consideration. As was stressed above, it is well- 
known that, due to the strong interference, the resonance reveals itself not only as a bump in the 
spectrum but also as a dip or a shoulder; likewise, the resonances appear in the OO"*"*" wave. Second, 
the interference effects are the source of ambiguities. It is worth noting that ambiguities in 00++ wave 
were successfully eliminated in the analyses of Refs. pQ, M, that was due to the simultaneous fit of 
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different meson spectra only. For the wave ^0^, the available data are not copious, hence one may 

suspect that the solution found in Ref. is not unique. 

The _fC-matrix re-analysis of the Ktt S-wave has been carried out in Ref. with a purpose: 
(i) to restore the masses and couling constants of the bare states for the wave ^O"*", in order to 

establish the gg-classification; 

{a) to find out all possible iiT-matrix solutions for the Ktt S'-wave in the mass region up to 2000 

MeV. 

The S'-wave Kt: scattering amplitude extracted from the reaction K^p — )■ K^n^n at small mo- 
mentum transfers is a sum of two components, with isotopic spins i and |: 

+ iAf =|As|e'*^ (3.17) 

where | As \ and (ps are measurable quantities entering the S'-wave amplitude The part of S- 

wave amplitude with the isotopic spin 1 — 3/2 has a non-resonance behaviour at the energies under 
consideration, so it can be parametrised as follows: 



1 - iPKTT (5)03/2(5) 



^ ^ ^ :^ ;7 , (s.is) 



where 03/2(5) is a smooth function and pK-nis) is the Kt: phase space factor. 

For the description of the A^J"^ amplitude, in |^ the 3x3 if-matrix was used, with the following 
channel notations: 1 = Ktt, 2 = Krj' , 3 = Ktttttt + multimeson states. The account of the channel 
Kf] does not influence the data description, since the transition Ktt — > Kr] is suppressed |3^, that 
also agrees with quark combinatoric results, see Table 3. In Ref. the fitting to the wave ^O"*" 
was performed with the parametrization of Kab given in (3.1). The analysed data for the reaction 
K~p K~-K^n were extracted with small momentum transfers {\t\ < 0.2 GeV^), and at the first 
stage the data were fitted to the unitary amplitude (3.1). At the next stage, the ^-dependence was 
introduced into the if-matrix amplitude. The amplitude Kn^t) Kn (7r(i) stands for virtual pion) 
is equal to A]/"^ = Kia{t) {I — ipK)~l: the parametrization of the matrix Kia{t) is written in (3.3). 

Coupling constants are determined by the rules of quark combinatorics, they are presented in Table 
3. In Ref. only the leading terms in the 1/A^ expansion were taken into consideration: in this 
case, coupling constants were fixed by the fit of the 00++ and 10++ waves, because is a common 
parameter for all the nonet members. 

The description of the ■iO+-wave has been performed under two assumptions, namely, with two- 
and three-pole structure of the wave in the mass region below 2000 MeV. 

In Ref. |3^ two solutions, A and B, were found for the wave ^0+; they differ at Mt^k > 1800 MeV 
only. Correspondingly, in Ref. the two two-pole if-matrix solutions, (A-1) and (B-1), have been 
obtained. The positions of amplitude poles are practically the same for both solutions; they are shown 
in Eq. (1.16). The description of data is shown in Fig. 12. The mass of the first resonance, see (1.16), 
does not differ strongly from that obtained in |3^, see (3.15), but the width of the resonance found 
in the i^T-matrix solution is twice as less. This follows from the doubling of poles in the ii'-matrix 
solution due to the correct account of the Kr/' threshold. The mass of the second resonance decreased 
for the iiT-matrix solution, as compared to the result of Ref. js^, in more than 100 MeV. 

The masses of bare kaon states related to the two-pole solution are given in Eq. (1.17). The mass 
of lightest state is 1200^^5o MeV, i.e. this scalar kaon is located in the same mass region as the other 
scalars, candidates for basic I^Pq nonet members. 

The description of data in the three-pole _fC-matrix fit is shown in Fig. 13. The region of high 
masses, Mktt > 1700 MeV is described by two poles in Solutions (A-2) and (B-2). However, the 
two-pole structure of the amplitude at Mktt > 1700 MeV did not influence the characteristics of two 
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low-located resonances — they are identical in Solutions (A-2) and (B-2). In Solution (B-3) the region 
Mk-k < 1600 MeV is described by two poles. Positions of bare states in Solution (B-3) are given in 
(1.18), the corresponding positions of poles are given in Eq. (1.19). 

Solid curves in Figs. 12 and 13 represent the description of the Kn wave by the unitary amplitude, 
and dashed lines represent the fits, where the t-dependence of the Ktt amplitude is taken into account. 
It is seen that the i-dependence allows us to get better description of phase shifts around 1700 MeV. 
It should be noted that in this region as well as in the mass region above 2000 MeV for Solution 
A, the data under consideration violate the unitary limit. It is hardly possible that rather strong 
violation of unitarity is a consequence of the amplitude i-dependence; it is more likely related to the 
underestimation of systematic errors in the partial- wave analysis of Ref. ]3^ in those regions. The 
inclusion of the f-dependence into fitting procedure does not affect strongly the masses of bare states 
and positions of pole of the amplitude. As a rule, the masses of bare states found in the i-dependent 
fit are less in 20-30 MeV than masses obtained in the t-independent fits. 
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Tabic 3 

Coupling constants for the transitions Kq two mesons and two mesons in the leading and 

next-to-leading terms of the 1/N expansion. 



Channel 


Couplings for 
leading terms 


Couplings for 
next-to-leading terms 


K+TT- 
















(cos e/V2 - VA sin 6)5^/2 


(V2cose - VAsine) g^^/2 




(sine/\/2-|-\/Acose) 


(V2sine - \/Acose) g^^/2 


K-K'' 







n~r] 


cose/v^ 


(\/2cose - \/Asine) g^^/2 


Tr~r]' 


sine/V2 


(V2 sin e - VA cos g^^/2 
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Fig. 12. Description of data in Ref. |^ in the two-pole Jf-matrix fit: Solutions (A-1) and (B-1). 
Solid curves correspond to the solution found for the unitary amplitude, dashed line stands for the fit 
with the i-dependent if-matrix. 
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Fig. 13. Description of data in Rcf. in the three-pole if-matrix fit: Solutions (A-2), (B-2) 
and (B-3). 

4 Propagator matrix: analysis of the {IJ^^ = 00^^)-wave 

Here we are summing up the results of the analysis of the 00++ 5- wave performed in ^, in 
terms of propagator matrix (D- matrix). The £)- matrix technique is based on the dispersion relation 
N / D-method, it allows us to reconstruct the amplitude, which is analytical on the whole complex s- 
plane. We discuss the effects which are due to the resonance overlapping and mixing: mass shifts and 
accumulation of widths by one of neighbouring resonances. The expansion of the physically observed 
states in a series with respect to initial (non-mixed) ones is performed. 

The investigation is done for the 00+"'" wave, though the method can be easily generalised for other 
waves using the technique developped in Refs. |4q| . 
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4.1 The mixing of two unstable states 

In the case of two resonances, the propagator of the state 1 is determined by the diagrams of Fig. 
14a. With all these processes taken into account, the propagator of the state 1 is equal to: 



Du{s) 



mi 



Bu{s) 



Bi2{s)B2l{s) 

ml - s - B22 (s) 



(4.1) 



Here mi and m2 arc masses of the input states 1 and 2, and the loop diagrams -By(s) are defined by 
Eq. (2.21), with the replacement g'^{s) gi{s)gj{s). It is helpful to introduce the propagator matrix 
Dij, where the nondiagonal elements -D12 = -D21 correspond to the transitions 1^2 and 2^1 (see 
Fig. 14b). The matrix reads: 

D 



Dn D12 
D21 D22 



{Mf - s){Mi - s) - B,2B2i 



21, 



-B12 
Ml-s 



Here the following notation is used: 



Bii{s) 



1,2. 



(4.2) 



(4.3) 



Zeros of the denominator of the propagator matrix (4.2) define the complex resonance masses after 
the mixing: 

n(s) = {Ml - s)iMl -s)~ B12B21 = . (4.4) 

Let us denote the complex masses of mixed states by Ma and Mb- 

Consider a simple model, where the s-dependence of the function -Bij(s) near the points s ^ M\ 
and s Mg is assumed to be negligible. Let Mf and B12 be constants, then one has: 



MXb = \{Ml + Mi) ± J i(M2 - M|)2 + 512^21 



(4.5) 



In the case, when the widths of initial resonances 1 and 2 are small (hence the imaginary part of the 
transition diagram B12 is also small), the equation (4.5) turns into the standard formula of quantum 
mechanics for the split of mixing levels, which become repulsive as a result of the mixing. Then 



D = 



COS 



- cos O sin i 
M'i-s 



Ml-s 



+ 



B 

sin ^ COS 6 
M'i-s 



- cos e/ sin t 

sin^ e 
Mi-s 



+ 



sm U cos t 



Mi 



I cos 



Ml-s 



, COS = - + - 



^\{Ml-MlY+B^2B2i 

(4.6) 



The states \A > and \B > are superpositions of initial levels, |1 > and |2 >, as follows: 
|A >= cos6'|l > -sin6'|2 > , IS >= sin6»|l > + cos6'|2 > . 



(4.7) 



In general, representation of states \A > and |B > as superpositions of initial states is valid, when 
one cannot neglect the s-dependence of functions Bij{s) and their imaginary parts are not small. 
Consider the propagator matrix near s = M\: 



D 



1 



n(s) 



M|(s) 



Sl2(s) 



B2i{s) Ml{s)-s 



-1 



W{Ml){M\-s) 



Mi{Ml)-M\ i?i2(Ml) 
B2i{Ml) Mf{Ml)-Ml 



(4.8) 

In the left-hand side of Eq. (4.8), the singular (pole) terms are the only surviving. The matrix 
determinant in the right-hand side of (4.8) is equal to zero: 



[Ml{M\) - Ml][Ml{M\) - Ml] - B,2{Ml)B2,{M\) = 



(4.9) 
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This equality follows from Eq. (4.4), which fixes n(Af^) 
mixing angle: 



0. It allows us to introduce the complex 



|A >= cos6'^|l > -sin6lA|2 > . (4.10) 
The right-hand side of Eq. (4.8) can be rewritten with the use of the mixing angle Oa, as follows: 

Na 



D 



-Ml 



Ml 



cos Ua 
■ sin 9a cos 9a 



cos Va sm t)A 
s\v? 9a 



(4.11) 



where 



Ml 



Ml 



Ml 



= ir(Ml) [^^^ - - ^^1' = 2Ml-M?-Mr = 2Ml-Mf-Mi ' ^''''^ 



We remind that in the formula (4.12) the functions Mj^(s), M|(s) and i?i2(s) are fixed in the point 

s = Ml. In the case under consideration, when the angle 9a is a complex quantity, the values cos^ 9a 

,2 , 



and sin 9a do not determine the probability of states |1 > and |2 > in \A >; indeed, the values 
\/Na cos 9a and —^Na sin 9a are the transition amplitudes \A >— > |1 > and \A >^ |2 >. Therefore, 
the corresponding probabilities are equal to | cos9a\'^ and | sin^^P- 

In order to analyse the content of the state \B >, an analogous expansion of the propagator matrix 
should be done near the point s = M%. After introducing 



\B >= sin 6*511 > +cos6Ib|2 > , 
we have the following expression for D in the vicinity of the second pole s 

Nb 



Ml: 



D 



Ml 



sin^ 9 b 



cos Vb sm Vb 
sin cos cos^ 6b 



(4.13) 



(4.14) 



where 
Nb 



II' {M 



[2M| 



Ml - Mi] , cos^ 9b 



Ml 



Ml 



2Ml - Ml - M| ' 



sin^ 9 b 



Ml 



Mi 



2M% - Ml - Mi ■ 



(4.15) 

In Eq. (4.15) the functions (s), Af|(.s) and Bi2{s) are fixed in the point s = Ml- 

If Bi2 depends weakly on s and one can neglect this dependence, the angles 9a and 9b coincide. 

But in general they are different. So the formulae for the propagator matrix differ from the standard 

approach of quantum mechanics by this very point. 

Another distinction is related to the type of the level shift afforded by mixing, namely, in quantum 

mechanics the levels "repulse" each other from the mean value l/2(£'i + E2) (see also Eq. (4.5)). 

Generally, the equation (4.4) can cause both the "repulsion" of masses squared from the mean value, 

l/2{Ml + Mi), and the "attraction". 

The scattering amplitude in the one-channel case is defined by the following expression: 



A{s) = gi{s)Dij{s)gj{s) . 
In the multichannel case, Bij (s) is a sum of loop diagrams: 

b,,(5) = ^b1;)(5), 



(4.16) 



(4.17) 



B^J^ is a loop diagram in the channel n with vertex functions g^^\ 
Partial scattering amplitude in the channel n is equal to: 



Qj"'^ and phase space factor p„. 



An{s)=gi-\s)D,j{s)g)">{s). 



in), 



(4.18) 
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4.2 The overlapping of a large number of resonances: construction of prop- 
agator matrix 

Consider the construction of the propagator matrix D for arbitrary number of resonances. The matrix 

elements, D^, describe the transition from the input state i (with the bare propagator (m| — s)~^) 
to the state j. They obey the system of linear equations as follows: 



Aj = DikBkj{s){m] - s) ^ + kjim] - s) ^ , 



(4.19) 



where Bij{s) is the loop diagram for the transition i j and 5ij is the Kronecker symbol. Let us 
introduce the diagonal propagator matrix d for input states : 

d = diag {{ml — s)~^, {m\ — s)~^, {m\ — s)~^ ' ' ') • 

Then the system of linear equations (4.19) can be rewritten in the matrix form as follows: 

b = DBd + d . 



One obtains: 

The matrix is diagonal, hence 

Ml 

b-' = 



D = — — . 

- {d^^ — 13) is of the form: 

- S -Si2(s) -Si3(s) 

B2i{s) Mi-s -B23{S) 



(4.20) 

(4.21) 
(4.22) 



Bsi{s) -Bs2{s) Ml-s 



(4.23) 



where Mf is defined by Eq. (4.3). Inversing this matrix, we obtain a full set of elements Dij{s): 



Dij{s) 



nW(s) 



(4.24) 



Here II^^) (s) is the determinant of the matrix D ^ , and Hj- ^' {s) is a matrix supplement to the 



element [D~^]ji, i.e. the matrix with the excluded j-th line and i-th column. 

The zeros of n(^)(,s) define the poles of the propagator matrix which correspond to physical 
resonances formed by the mixing. We denote the complex resonance masses as: 



s = Ml, 



Ml, 



Ml,. 



(4.25) 



Near the point s = M\, one can leave in the propagator matrix the leading pole term only. This 
means that the free term in Eq. (4.21) can be neglected, so we get a system of homogeneous equations: 



Dik{s) fj-i -B] = 

V / kj 



(4.26) 



The solution of this system is defined up to the normalization factor, and it does not depend on the 
initial index i. Then the elements of the propagator matrix may be written in a factorized form as 
follows: 



Na 



Ml 



Oil, 



a\(X2, oi\az, 



(4.27) 
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where Na is the normahzation factor chosen to satisfy the condition: 

al + al + al + ... + alf = 1. (4.28) 

The constants ai arc the normahzcd amphtudcs for the transitions: resonance A — > state i. The 
probabihty to find the state i in the physical resonance A is equal to: 

Wi = |Qi|2 . (4.29) 

Analogous representation of the propagator matrix can be also done in the vicinity of other poles: 

DiP{s^Ml)=NBj^^, Df\s^Ml) = Nc^p^-^ •••• (4.30) 

Coupling constants satisfy normalization conditions similar to that of Eq. (4.28): 

+ . ..+Pl = l, 7?+7l + ---+7^ = l, •••• (4.31) 

In general case, however, there is no completncss condition for the inverse expansion: 

a' +/?'+7' + ---?^l- (4.32) 

For two resonances, it means that cos^ Qa + sin^ 0b ^ 1. Still, let us remind that the equality in the 
inverse expansion, which is relevant to the completness condition, appears in the models, where the 
s-dependence of loop diagrams is neglected, see Eqs. (4.5)-(4.7). 

4.3 Full resonance overlapping: the accumulation of widths of neighbour- 
ing resonances by one of them 

Let us consider two examples which describe the idealized situation of a full overlapping of two or 
three resonances. In these examples, the effect of accumulation of widths of neighbouring resonances 
by one of them can be seen in its original untouched form. 

a) Full overlapping of two resonances. 

For the simplicity sake, let Bij be a weak s-dependent function, and Eq. (4.5) may be used. We 
define: 

Ml = Ml- iMnTi , M| = M| - zMflFa , (4.33) 

It is possible that ReBi2{M'^) can be equal to zero at positive gi and 92, if the contribution from the 
integration region s' < cancels the contribution from the region s' > M^. In this case 

Si2(M2) ^ ig^{Ml)g2{MDp{Ml) = iMR^/f^2. (4.35) 

Substituting Eqs. (4.33)-(4.35) into Eq. (4.5), one has: 

M\^Ml- iMniVi + T2) Ml^Ml. (4.36) 

Therefore, after mixing, one of the states accumulates the widths of primary resonances, — » F1+F2, 
and another state becomes quasi-stable particle, with Fb — > 0. 

b) Full overlapping of three resonances. 
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Consider the equation 

n(3)(s) = (4.37) 
at the same approximation as in the above example. Correspondingly, we put: 

ReBabiMl) ^ , (a 7^ &); Mf = Ml-s- iMnVi = x - i^i . (4.38) 

A new variable, x = — s, is used, and we denote MaTi = ji. Taking account of BijBji = —^i'yj 
and B12B23B31 = —1717273, we can re-write the equation (4.37) as follows: 

x^ + x^ (7:71 + 272 + 473) = . (4.39) 

Therefore, at full resonance overlapping, one obtains: 

Ml^Ml,-iMn{T,+T2+Ts), ^ M% , ^ M^. (4.40) 

The resonance A has accumulated the widths of three primary resonances, and the states B and C 
became quasi-stable and degenerate. 



4.4 The resonances /o(1300), /o(1500), /o(1530lio) and /o(1780) 

The ii'-matrix analysis provides a basis for the investigation of mixing phenomenon in the scalar 
sector. Propagator matrix technique used at the next stage of the analysis allows us to restore 
correctly the contribution from the real parts of loop diagrams, Bij (s) , thus having correctly calculated 
contributions of input states to the formation of physical resonances. 

The resonance mixing in the region 1200-1600 MeV may be considered in the two-channel approx- 
imation, for the quark-hadron duality justifies the analysis with the use of quark channels nh and ss. 
Correspondingly, 

Bij{s) = cosipi cos ipjBl""'\s) + sinipi sin<^jB|**^(s) , (4-41) 

where i, j run over 1,2,3,4, with the following notations for states: 1 = l^Po{nn rich), 2 = 
2^Pa{nn rich), 3 = gluonium and 4 = 2^Pq[ss rich). The quark states are usually described by 
the light cone variables. Then, 

,(nn). ^ 1_ dx [ ^ 9i{s')gj{s') , 



(.) -J^liJ "'•'^'-'^^l'' - 4™') . (4.42) 

Here s' = {m? + /c^)/x(l — x) and m is the mass of nonstrange quark. The factor 2(s' — 4m^) appears 
due to the presence of the quark spin variables: Tr[{k + m){—p + k + m)] — 2{s' — 4rn^). Analogous 
expression, with the replacement m — > m^, determines bI^^\s). 



The simplest parametrization of vertex function for the transition input state i quarks is: 



9i{s) = 7i 



fc2 -h CT, k'^ + ai + h 



(4.43) 



Here fc^ = s/4 — m? and fc^ = to^/4 — m^, where m is the constituent quark mass equal to 350 MeV 
for nonstrange quark and 500 MeV for strange one, and rUa is the mass of input state. 

For the first state, l^Po{nn rich), and for the gluonium we put ^1=^3 = 0. The second state is 
the radial excitation, 2^Pci{n.h rich), and its wave function is orthogonal to the groimd state. This 
means that real part of the function Bi2{s) must tend to zero at s close to resonance masses. Such an 
orthogonalization has been performed at the point ^/s = 1.5 GeV, thus determining the value of the 
coefiicient d2. Vertex functions for the same nonet members are equal to each other, so 52(5) = 54(5)- 
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Parameters ttIq, 7^, h and tXa are defined by the masses and widths of physical resonances. However, 
the mass iJia can be approximately fixed by the ii'-matrix pole: /i^ ~ — Re BaailJ^)- It should be 
underlined that is the mass of a pure gluonium which is a subject of the Lattice QCD. 

Positions of amplitude poles and the masses of input states found in Ref. by fitting the 00++ 
wave are shown in Table 4. Relative weight of the primary state in the physical resonance A is defined 
by Eq. (4.29): calculated in such a way the probabilities Wa for the resonances under investigation 
are shown in Table 4. 

As was stressed above, in order to make comparisons with the QCD calculations, one should 
separate the contributions from large and small distances, that is, to take into account the short-range 
interaction component, r < ro ~ Rconfinement, and eliminate the contribution from large r. Therefore, 
in the calculation of masses which might be compared with the results of the QCD-motivated models, 
we should make a replacement in the amplitude of the 00++ wave as follows: 

00 

B^bis) ReB^bis.kl) = P I ^ 9A^M^:M^ 2is' - Arn^) . (4.44) 

J TT S' — S 

The poles of the amplitude re-determined in this way provide the masses which are related to the 
interaction at r < 1/fco- One must compare with the quark model results the masses obtained with a 
cutting, of the order of ^ 0.125 (GeV/c)^, that corersponds to the account of the quark interaction 

•^t r ^ 1 flTL ^ Rconfinement • 

For Solution I we get (the values are given in GeV): 

l^PQ{ssrich) l^PQ{nnrich) 2^Pf){nn rich) 2'^ Pq{ss rich) 

m{k^ =0)^ ^l^'-^' 0.720 1.360 1.577 1.791 (4.45) 

m(fcg = 0.125) 0.730 1.340 1.560 1.780 

m{kl -> 00) = rua - 1.457 1.536 1.750 

In Solution II: 



l^PQ{ssrich) l^Po{nnrich) 2^Po{nnrich) 2^PQ{ssrich) 

m(fc2 = 0) /i^^'^'''^ 0.720 1.357 1.585 1.734 

m(fc2= 0.125) 0.735 1.340 1.570 1.725 

m(fc^ -^00) ^rua - 1.107 1.566 1.702 



(4.46) 



The lightest gg-state, /q°'''^(720), has not been included into the mixing machinary in Refs. 
In Eqs. (4.45) and (4.46) the mass corrections for this state have been evaluated as ma(fcQ) ~ 

— Re BaainT-aj ^o)- This approximate equality is due to a comparative smallness of the nondiagonal 
loop diagrams. 

Equations (4.45) and (4.46) prove that the values uiaikQ — 0.125 GeV^) slightly differ from ^^'""'^) 
while the differences from the input masses ma can be significant. This means that the /^-matrix 
analysis provides approximately correct meson characteristics which may be compared with the quark 
model results. On the contrary, one should compare with the results of Lattice QCD the values of 
input masses rUa, which can noticeably differ both from the masses of bare states, n^""^^, and from 
those of real resonances. 

4.5 Dynamics of glueball mixing with gg-states 

To find out the glueball mixing with qq- states, let us make a replacement in loop diagrams of the 
propagator matrix: 

ga{s) ^ ^9a{s) , (4.47) 
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with the factor ^ changing in the interval 0<^<1- At^^O the mixing is switcliecl off, and 
the amphtude has the poles at s ~ m^. Figure 15 demonstrates the position of poles at different ^ 
for Solutions 1 and II. With increasing ^, the poles move from the real axis to the lower part of the 
complex plane. Let us discuss in detail the dynamics of the pole movement for Solution II. 

At ^ = 0.1 — 0.5 the glueball state of Solution II is mainly mixed with the state 2^ PQ{nn rich) , 
while at ^ = 0.8 — 1.0 the mixing with the state l^Pf){nnrich becomes important. As a result, the 
state, which is the glueball descendant, sunk rather deeply into complex plane, acquiring the mass 
M = 1450 — i450 MeV, and the glueball component of this broad resonance is 47%. Likewise, in 
Solution I the broad resonance is the gluon descendant as well. 

The hypothesis about strong mixing of the gluonium with the qq-states was raised formerly. But 
the attempts to restore quantitatively the picture of mixing within the standard quantum mechanics 
approach failed, for two phenomena have been lost: 

(1) The qq/ glueball mixing described by the D-matrix may lead not only to the repulsion of levels, 
that follows from the standard quantum mechanical approach, but to the attraction of levels as well. 
The latest effect is caused by the presence of imaginary parts of the loop diagrams Bab, and it is 
important that Im Bab is not small near 1500 MeV. 

(2) The resonance overlapping leads to the repulsion of poles located on the imaginary axis of 
masses, and one resonance accumulates the widths of the others. 

This very type of the mixing occurred at 1500 MeV, and a large width of one of the resonances is 
its inevitable consequence. 

It is also natural that the broad resonance itself is the gluonium descendant, for the gluonuim 
mixes without any significant suppression with the nearby gq-states, both of them being dominantly 
nonstrange ones. 
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Tabic 4 



Masses and mixing angles of the input states, 
the content of physical states and positions of poles of the 00++-amplitude (masses in 



Solution I 




l'^PQ{nn rich) 


2'^PQ{nn rich) 


Gluonium 


2''^ Pq{ss rich) 




01 = 18° 


02 = -6° 


03 = 25° 


04 = 84° 




mi = 1.457 


m2 = 1.536 


ms = 1.230 


m4 = 1.750 


iy[/o(1300)] 


32% 


12% 


55% 


1% 


1.300-^0.115 










Vt^[/o(1500)] 


25% 


70% 


3% 


2% 


1.500 -i0.065 










VK[/o(1530)] 


44% 


24% 


27% 


4% 


1.450 - i0.450 










W^[/o(1780)] 


1% 


1% 


0% 


98% 


1.780 - i0.085 










Solution II 




l'^Po{nn rich) 


2^Po{nn rich) 


Gluonium 


2^Po(ss rich) 




01 = 18° 


02 = 35° 


03 - 25° 


(l)A = -55° 




mi = 1.107 


m2 = 1.566 


ms = 1.633 


m4 = 1.702 


T^[/o(1300)] 


35% 


26% 


38% 


0.4% 


1.300 -i0.115 










VF[/o(1500)] 


1% 


64% 


35% 


0.4% 


1.500 -i0.065 










W^[/o(1530)] 


12% 


41% 


47% 


0.3% 


1.450 -i0.450 










W^[/o(1780)] 


0.1% 


0.2% 


0.2% 


99.5% 


1.750 - iO.lOO 
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Fig. 14. Diagrams describing the propagation functions Du (a) and £'12 (b) for the interaction 
of two bare states. 
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M (GeV) 



Fig. 15. Variation of the resonance pole positions plotted on the complex-M plane at different 
coupling constants ga £.ga- 

5 Conclusion 

Deconfinement of quarks from the excited gg-levels is going in two stages: 

1) Unavoidable production of qq-pairs, which form two or more white states (hadrons). 

2) The flying away of the produced hadrons, their interaction and, as a result, the mixing of neighbour- 
ing gg-levels that leads to the formation of a broad locking state, which plays the role of a dynamical 
barrier for neighbouring levels. 

That is the /^-matrix analysis together with the dispersion N/D method, which are summoned to 
decypher the second stage of the deconfinement. Analysis of 00++-wave carried out in the /iT-matrix 
or propagator matrix techniques demonstrated that the lightest scalar glueball, being near the states 
l^Poqq and 2'^Poqq, turned after the mixing into a broad state, with r/2 ~ 500 MeV. This broad state 
/o(1530^25o) carries about a half of the scalar gluonium component, while the other components of 
the broad resonance are l^Poqq and i^P^qq. 

It looks like the waves 00 ^ and 02++ behave similarly |5^, that allows us to believe that physics 
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of highly excited states is tightly connected with the study of broad states, being an important and 
inevitable step in the search for exotic hadrons and the investigation of confinement phenomenon. 



Thanks are due to L.G. Dakhno and V.A. Nikonov for help. This investigation is supported by RFBR 
grant N 96-02-17934 and the INTAS-RFBR grant N 95-0267. 
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